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ABSTRACT 

The driving mechanism of jets and outflows in star formation process is 
studied using three-dimensional resistive MHD nested grid simulations. Start- 
ing with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, 
we calculated cloud evolution from the molecular cloud core (n c = 10 4 cm -3 , 
r = 4.6 x 10 4 AU) to the stellar core (n c = 10 22 cm 3 , r ~ 1-R©), where n c and r 
denote the central density and radius of each object, respectively. In the collaps- 
ing cloud core, we found two distinct flows: Low- velocity flows (~5kms _1 ) with 
a wide opening angle, driven from the adiabatic core when the central density 
exceeds n c > 10 12 cm~ 3 , and high- velocity flows (~30kms _1 ) with good collima- 
tion, driven from the protostar when the central density exceeds n c > 10 21 cm" 3 . 
High-velocity flows are enclosed by low-velocity flows after protostar formation. 
The difference in the degree of collimation between the two flows is caused by 
the strength of the magnetic field and configuration of the magnetic field lines. 
The magnetic field around an adiabatic core is strong and has an hourglass con- 
figuration; therefore, flows from the adiabatic core are driven mainly by the 
magnetocentrifugal mechanism and guided by the hourglass-like field lines. In 
contrast, the magnetic field around the protostar is weak and has a straight con- 
figuration owing to Ohmic dissipation in the high-density gas region. Therefore, 
flows from the protostar are driven mainly by the magnetic pressure gradient 
force and guided by straight field lines. Differing depth of the gravitational po- 
tential between the adiabatic core and the protostar cause the difference of the 
flow speed. Low-velocity flows correspond to the observed molecular outflows, 
while high-velocity flows correspond to the observed optical jets. We suggest that 
the outflow and the jet are driven by different cores, rather than that the outflow 
being entrained by the jet. 
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Introduction 



Since Snell discovered an outflow from a protostar in 1980 (ISnell et al. Ill980l ). over 300 
outflows have been found in many star-fo rming regions rtWu et al. Il2004h. Recently, ou tflows 
from brown dwarf (jWhelan et al. 1120051 1 and O-type stars (IBeuther fc Shepherd! 120051 ) were 
discovered. These observations indicate that outflow is ubiquitous in the star formation 
proce ss. Observations a l so show that outflows h ave two or more distinct velocity components 



le.g. 



Hirth et al. 1119971 ; iPyo et al. II20031 . 120051 ) . Typically, an outflow is composed of a low- 
velocity com ponent (LVC) wit h 10-50 km s _1 and a high- velocity component (HVC) with 
~100kms _1 (IPyo et al. 1 120051 ) . Outflows have a collima tion factor (i.e., l ength/width or 
major/minor radius of the outflow) of H/R = 3 to > 20 (jArce et al. 1120061 ). where H and 

R are the height and radius of the outflow , respectivel y. There is a clear trend toward 

higher collimation at higher outflow velocity (jArce et al. 1120061 ) . Since outflows have various 
morphological and kinematical properties, they cannot be explained by a single-class model. 
The flows that originated from the protostar are typically classified into two types: molecular 
outflow (hereafter, outflow) obse rved mainly with CO molecules, and o ptical jet (hereafter 
jet) observed by optical emission (jArce et al. 1120061 ; iPudritz et al. 1120061). Outflow s observed 



by CO line emission exh ibit a wide opening angle (e.g., iBelloche et al. 1 120021 ) and slow 



velocity (10 — 50kms ; iRicher et al. I I200CH : lArce et al. 1120061). Jets observed by optical 



emission, however, exhibit good coll imation (e.g., iHirano et al. 1 120061 ) and high velocity 



100 — 500 km s x ; lArce et al. 1 120061 ) . Observations indicate that around each protostar, 



high-speed jets with a n arrow opening angle a re enclosed by a low-velocity outflow with a 



wide opening angle (e.g.. iMundt fc Friedlll983l ) 



What mechanism drives the outflow and the jet in the star formation process? Dur- 
ing star formation, neither radiation nor thermal pressure from the protostar can supply 



suffic i ent kinetic energy a nd momentum for driving the jet and the outflow (jWu et al. 



2004 : IPudritz et al. II2006I ). These flows (outflows and jets) are considered to be driven 



by release of the gravita t ional energy med iated by the Lorentz and centrifugal forces (cf. 
Blandford fc Paynelll982l ; IShu et al. Ill994l ). However, we have not understood what drives 
these flows during star formation. It is difficult to directly observe the driving points of the 
outflow and the jet, because they are embedded in the dense cloud at several AU from the 
protostar. Furthermore, we cannot analytically investigate the driving mechanisms of the 



- 3- 



jet and the outflow in detail, because they have a large spatial scale and density contrast, 
and these flows are influenced by closely related Lorentz, centrifugal, thermal pressure gra- 
dient, and gravitational forces. Therefore, numerical simulation is the most effective tool for 
investigating the driving mechanisms of the jet and the outflow. 



Blandford fc Payne! (119821 ) showed that centrifugally driven outflow from the disk is 



possible if the poloidal comp onent of the magne t ic fiel d mak es an angle of less than 60° 
(i.e., the disk wind model). iPudritz fc Norman! (119831 . Il986l ) proposed that this models 
as the mechanism for the protostellar jet. Subsequently many numerical simulations (e.g., 
Kudoh &: Shibata Hl997aU bl). which showed t hat disk could driye the j et confirmed the disk - 



wind mech a nism around the protostar (see, iKonigl fc Pudrita |2000| ; IPudritz et al. 1120061 ). 



Shu et al. I (119941 ) proposed the X-wind model, in which the jet is driven in close proximity 
to the protostar. Both the disk wind and the X-wind are driven magnetocenrifugally from 
open magnetic field lines anchored on rapidly rotating circumstellar disks. The difference 
between these models lies in where the field lines are anchored: near the radius of magne- 
tosph erical truncation on the disk for X-wind and over a wider range of disk radii for disk 
wind ( jShang et al. II2006I ). The jet speeds derived from both X-wind and disk wind theories 
correspond to the Kepler speed at their driving point on the circumstellar disk. Thus, these 
models cannot explain the two velocity components of jet and outflow. Many people believe 
that outflow (i.e., slow speed and wide opening angle flow) is entrained by the jet (i.e., high 
speed and narrow o pening angle flow) , and that the jet is driven by either the disk wind 



or the X-wind (e.g.. lArce et al. 1 120061 ) . However, almost all simulations, except for a few 
studies, calculate the evolution of the circumstellar disk (i.e., jet driving) after the protostar 
grows up sufficiently, and provide the idealized boundary, density profile of the disk, distri- 
bution of the magnetic field and velocity, accretion rate, etc. After the protostar is formed, 
we cannot judge whether these configurations are realized. Therefore to acquire adequate 
information about the protostar and to understand the driving mechanism of the jet and 
the outflows, we should calculate cloud evolution from the molecular cloud that forms the 
protostar. 

Simulations of cloud evolution from the molecular cloud core to stellar core forma- 
tion show that the low-velocity flow is driven from the adiabatic core formed before pro- 
tostar formation, and that the high-velocity flow is driven from another core. Assuming 
spherical symmetry, the evol ution from molecular cloud to stellar co re has been investi- 
gated by many authors (e.g., Larson 1969 ; Winkler fc Newman 1980al lbl; iMasunaga et al. 



19981 ; IMasunaga fc Inutsuka 1 120001 ). and they determined the density, velocity, and tem- 
perature structures of the collapsing cloud and stellar core. They have shown that the 
adiabatic core (or the first core) with a radius of ~1 AU forms when the central density 
reaches n c ~ 10 11 cm -3 , because the gas becomes optically thick and the thermal pressure 
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increases. iBatd ( 119981 ) and IWhitehouse fc Batd (120061 ) calculated the stellar core formation 
from the molecular cloud core using SPH simulations. However, the outflow and the jet do 
not appear in their calculations because they ignored the magnetic effect. The evolution 
of the magnetized cl oud until the fo r mation of th e a diabat ic core (hereafter, the fir s t core ) 



was investigated bv Machida et al. I (12004 l2005al ib . l2006af). iHosking & Whitworthl (J2004J) 



Matsumoto fc Tomisakal (120041 ). Izierierl (l2005h. and 



pears after t he first core formation i n lMachida et al. 



(120041) . and iFromang et al. I (I2006h . 



romang et a, 



M 



(2004, 2005b) 



Tomisakal Jl99sl . I2OO0I . l2002h . and baneriee k Pudritz 



( 20 06|). The outflow ap- 



Matsumoto &: Tomisaka 



( 120061 ) investigated the evolution of the magnetized clouds up to protostar formation. They 
calculated the cloud evolution from the molecular cloud (n c ~ 10 2 — 10 6 cm -3 ) to protostar 
formation (n c ~ 10 22 cm -3 ). In their calculations, two distinct flows appeared: low- velocity 
flow (v ~ 2kms _1 ) driven from the first core, and high-velocity flow (v ~ 30 kms -1 ) driven 
from the second core (i.e., the protostar). They expected that low- velocity flow would be 
observed as molecular outflow and high-velocity flow as an optical jet. They adopted ideal 
MHD approximation, which is valid in the low-density gas region ( n < 10 12 cm~ 3 ) ; how- 
ever, it is not valid in the high-density gas region (n > 1 '>'-'<■•" 



found sign ificant ma g netic flux loss for 10 cm < n < 10 cm by Ohmic dissipation. 



10 i2 cm" 3 ). 

15 , 



Nakano et al. I (120021 1 



Therefore, iTomisakal Jl998l . boool . I2OO2I ) and iBanerjee fc Pudritzl (120061 ) overestimated the 
magnetic flux of the cloud, particularly in the high-density gas region. 

There are two predictions concerning driving outflows and jets: one is that the outflow is 
entrained by the jet, and another is that the outflow and jet are driven from different objects. 
These predictions are mutually conflicting. In this paper, we calculate cloud evolution from 



10 4 cm 



4.6 x 10 4 AU) to stellar core formation 



10 22 cm 3 , r c ~ l-R©) using three-dimensional resistive MHD nested grid method, 



the molecular cloud core 

(n c 

study the formation process of jets and outflows, and show the driving mechanisms of these 
flows. The numerical method of our computations and our model framework are given in 
§2, and the numerical results are presented in §3. We discuss the driving mechanism of the 
outflow and the jet in §4, and summarize our conclusions in §5. 



2. Model 



Our initial settings are almost the same as those of iMachida et al. I J2006al Jbl. [2007l l To 
study cloud evolution, we use the three-dimensional resistive MHD nested grid code. We 
solve the resistive MHD equations including self-gravity: 



dp 
di 



V ■ (H = 0, 
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dv 

dB 



-VP 



1 

47T" 



B x (V x B) 



pV0, 



dt 



V x (v x B) +r/V 2 B, 
V 2 = AnGp, 



(2) 

(3) 
(4) 



where p, i>, P, B, 77, and denote the density, velocity, pressure, magnetic flux density, 
resistivity, and gravitational potential, r espectively. The res istivity (77) in equation (jHJ) is 
a function of density and temperatu r e (INakano et al. 112002 ). For resistivity (77), we use 
the value adop ted in iMachida et al. 1 (l2006bL boOTf ). To mimic the temperature evolution 
determined by iMasunaga fc Inutsuka I J2000h . we adopt the piece-wise polytropic equation 
of state as 



P 



C 2 S P 
4Pc 

C 2 sPc 

C 2 s Pc 



7/5 
7/5 
7/5 

-13 , 



1.1 



1.1 



5/3 



P < Pc, 
Pc< P < Pd, 
p d < p< p e , 

P > Pe, 

10 11 



(5) 



where c s = 190 ms , p c = 3.84 x 10~ i3 g cm~ 3 (n c = 10 il cm~ ii ), p d = 3.84 x 10~" g cm" 
(rid = 10 16 cm~ 3 ), and p e = 3.84 x 10~ 3 g cm -3 (n e = 10 21 cm~ 3 ). For convenience, we 
define "the protostar formation epoch" as that at which the central density (n c ) reaches 
n c = 10 21 cm -3 . We also call the period for which n c < 10 11 cm -3 "the isothermal phase," 



the period for which 10 cm < n c < 



10 16 cm" 3 



which 10 16 cm 3 < n r 



'the adiabatic phase," the period for 
< 10 21 cm -3 "the second collapse phase," and the period for which 



n c > 10 21 cm 



-3 » 



the protostellar phase.' 



In this paper, we adopt a spherical cloud with critical Bonnor-Ebert (]Ebertlll955l : iBonnor 
19561 ) density profile having p Cj o = 3.841 x 10 _20 g cm -3 (n C) o = 10 4 cm~ 3 ) of the central 



(number) density as the initial condition. The critical radius for a Bonnor-Ebert sphere 
R c = 6A5c a [4irGp B E(0)]~ 1/2 corresponds to R c = 4.58 x 10 4 AU for our settings. Initially, 
the cloud rotates rigidly (Qo) around the z-axis and has a uniform magnetic field (Bq = 17/iG) 
parallel to the z-axis (or rotation axis). To promote contraction, we increase the density by 
70% from the critical Bonnor-Ebert sphere. 

The initial model is characterized by a single one non-dimensional parameter u. This 
parameter is related to the cloud's rotation rate, and is defined using a central density p as 



Li) 



n /{47iGp Q ) 



1/2 



(6) 



The model parameters u>, magnetic field B , angular velocity Q Q , total mass inside the critical 
radius M , and the ratio of the thermal ao, rotational /3q, and magnetic 70 energies to the 
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gravitational energy^ are summarized in Table HJ 



We adopt the nested grid method (for details, see iMachida et al. I l2005al l2006al ) to 



obtain high spatial resolution near the center. Each level of a rectangular grid has the same 
number of cells (64 x 64 x 32), although the cell width h(l) depends on the grid level I. The 
cell width is reduced by a factor of 1/2 as the grid level increases by 1 (I — > l + l). We assume 
mirror symmetry with respect to z = 0. The highest level of a grid changes dynamically. We 
begin our calculations with four grid levels (I = 1, 2, 3, 4). The box size of the initial finest 
grid I = 4 is chosen to be 2R C , where R c denotes the radius of the critical Bonnor-Ebert 
sphere. The coarsest grid I = 1, therefore, has a box size of 2 4 R c . A boundary condition 
is imposed at r = 2 4 R c , such that the magnetic field and am bient gas rotate at an angular 
velocity of f2 (for details, see iMatsumoto fc Tomisakall2004l ). A new finer grid is generated 
whenever the minimum local Aj becomes smaller than 8 h(l max ). The maximum level of grids 
is restricted to Z max = 30. Since the de nsity is highest in the finest grid, the generation of a 
new grid ensures the Jeans condition of iTruelove et al. I (119971 ) with a margin of safet y facto r 
of two. We adopted the hyperbolic divergence B cleaning method of iDedner et al. I (120021 ) . 



3. Results 



Assuming spherical symmetry, many authors investigated the evolutio n from molec- 
ular cloud to stellar core using radiation hydrodynamic calculations (e.g., iLarsonl Il969 



Winkler fc NewmaJl980al lbl: lMasunaga et al. ll99alMasunaga fc Inutsuka koOoh ^We briefly 
summarize the evolution of the collapsing cloud core, according to their calculations. The 
molecular gas obeys the isothermal equation of state with a temperature of ~10K un- 
til n c ~ 5 x 10 10 cm~ 3 (isothermal phase), then the cloud collapses almost adiabatically 



cm 3 < n c < 10 16 cm 3 (adiabatic phase). A quasi-static core (i.e., the first 
10~ 2 M Q forms in the adiabatic phase. Subsequently dissociation of molecular 

16 , 



for 5 x 10 10 
core) with ~ 

hydrogen when the density exceeds n c > 10 16 cm -3 triggers further gravitational collapse 
(i.e., second collapse). Finally, the second core (i.e., protostar) with ~10 -3 M forms at 



t 10 cm . The protostar increases its mass for the subsequent gas accretion phase. 
In this study, we calculated the cloud evolution from the molecular cloud core wi th n, 



10 cm to form a protostar with n c ~ 10 cm" " 3 . IMachida et al. I (booi l2005al fbL l2006ah 
investigated cloud evolution in the isothermal and early adiabatic phases. The evolution of 
the magnetic field and angular velocity in the collapsing cloud for 10 4 cm -3 < n c < 10 23 cm -3 



1 Representing the thermal, rotational, magnetic, and gravitational energies as U, K, M, and W, the 
relative factors against the gravitational energy are defined as «o = po — K/\W\, and 70 = M/[W|. 
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was shown by iMachida et al. I (120071 ). Therefore, this paper focuses mainly on the cloud 
evolution after the first core formation (n c > 10 11 cm -3 ) because we are interested in the 
driving mechanism of the outflow and the jet. We parameterized angular velocity, magnetic 
field strength, and the ratio of thermal to gravitational energy of the initial cloud for 80 
models. Cloud evolutions are typically classified into three types when the outflow or jet 
appears in the collapsing cloud. We also found that cloud evolution depends on the ratio 
of the angular velocity to the magnetic field (Q /B ), but does not depend on the ratio of 
thermal to gravitational energy, especially after the first core formation (n c > 10 n cm~ 3 ). 
Observations indica te that molecu l ar cloud cor es have large magnetic energies and small 
rotational energies (jCrutcherlll999l ; ICasellil 120021 ) . In these types of clouds, c l oud evolutio n 
is sensitive to the initial angular velocity, as shown in IMachida et al. I (l2005al . l2006al . 120071 ). 
In the following, we show only three typical models in which the initial clouds have the same 
magnetic field strengths and ratios of thermal to gravitational energy, but different angular 
velocities (slow, moderate, and rapid rotation rate). 



3.1. Cloud Evolution with Slow Rotation 



3.1.1. First Core and Protostar Formation 



Figure [JJ shows the cloud evolution for model SR after the first core formation. Model 
SR has a parameter of uj = 0.003, which means that the initial cloud has a small rotation 
energy; the ratio of the rotational to gravitational energy is (3q = 3 x 10 -5 (see Table [TJ). 
Table [JJ shows that, in model SR, the rotational energy is much smaller than the thermal 
and magnetic energy in its initial state. 

The density (false colors) and velocity (arrows) distribution around the first core are 
plotted in the upper and middle panels of Figure [TJ The plasma beta and magnetic fields 
around the first core are plotted in the lower panels of Figure [TJ We plotted shocked regions 
as dotted lines in th e upper and middle panels in Figure [TJ Since the first core is surrounded 
by the shock layer (IMasunaga fc Inutsuka 1 12000| ) . dotted lines indicate the first core. For 



model SR, the first core forms when the central density reaches n c ~ 10 cm" . The first 
core at its formation epoch has a mass of 5.1 x 10~ 3 M Q and a radius of R ~ 0.71 AU, which 
is equivalent to or slig htly larger than that expect e d by one-dimensional radiatio n hydro- 
dynamic calculations (IMasunaga fc Inutsuka Il2000l ). iMatsumoto fc Hanawal (120031 ) showed 
that a larger first core forms in a rapidly rotating cloud because the centrifugal force sup- 
presses the cloud collapse and the shock occurs in the earlier adiabatic phase. However, the 
effect of centrifugal force at the first core formation epoch is small in model SR because a 
slowly rotating cloud is adopted as the initial state. The magnetic effect slightly increases 
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the mass and size of the first core, as shown in iMachida et al. I (j2005al ). 



Figures [TJa and b show the density and velocity distribution on z = and y = cut 
planes at n c = 6.5 x 10 16 cm -3 . These panels indicate that the first core has a nearly spherical 
shape. To evaluate the core shape, we define the oblateness as e Q b = (hih s ) l l 2 /h z , where hi, 
h s , and h z are the major axis, minor axis, and z-axis, respectively, derived from the moment 
of inertia for the high-density gas of p > 0.1p c according to Matsumoto & Hanawa (1999). 
Figure [2]a shows the evolution of the oblateness as a solid line. The oblateness increases as 
the cloud collapses in the isothermal phase (n c < 10 14 cm~ 3 ), and has a peak of e Q ^ ~ 10 
(i.e., the ratio of the radial to vertical scale is about 10) at n c ~ 2 x 10 8 cm -3 . Therefore, 
the central region has a disk-like structure at this epoch. Since the cloud rotates slowly, the 
disk is formed mainly by the Lorentz force. Note that the magnetic energy is comparable to 
both the thermal and gravitational energy at the initial state (see Tabled]). The oblateness 
oscillates around £ D b — 10 for 2 x 10 8 cm -3 < n c < 10 12 cm~ 3 , and then begins to decrease 
(Fig. [T]a). This decrease is caused by an increase in thermal energy. When the gas density 
reaches n c ~ 10 11 cm -3 , the cloud collapses adiabatically and thermal pressure increases. 
The oblateness becomes £ Q b — 1 at the first core formation epoch (n c ~ 10 14 cm~ 3 ), which 
indicates that the first core has a nearly spherical shape, as shown in Figure [T]a and b. 

Figure [T]c shows the plasma beta (/5 P = B 2 /Siic 2 p) around the first core on the y = 
cut plane. This panel shows that the magnetic energy inside the first core is extremely 
small (j3 p ~ 10 — 10 4 ), while that outside the first co re is comparable to the thermal energy 



(P p < 1). As shown in IMachida et al. I (l2006bl . 120071 ). the magnetic field is largely removed 



from the cloud core by Ohmic dissipation for 10 12 cm -3 < n c < 10 15 cm~ 3 . Figure [2]6 shows 
the evolution of the magnetic flux density at the center of the cloud (B c ) for model SR (solid 
line). The growth rate of the magnetic flux density becomes small for 10 12 cm~ 3 < n c < 
10 15 cm -3 . The solid line in Figure [2]c shows the magnetic flux density normalized by the 

1 /2 

square root of the gas density at the center of the cloud [B c j p c ) for model SR. In the ideal 

i li 

MHD regime, the evolution of B c / p c depends on the geometry of the collapse. This value, 

1 /2 

B c / p c , remains constant after a thin disk forms because the magnetic field increases in 
proportional to the s quare root of the density when the disk-like cloud collapses (B oc p 1 / 2 ; 
Scotto fc Black 1980 ). On the other hand, B C J pj 2 increases in pro portion to p 1 / 6 when the 
cloud collapses spherically (for details, see IMachida et al. Il2006al ). Thus, a rapid drop of 



B c /p\ /2 for 10 12 cm 3 < n c < 10 15 cm 3 in dicates the removal of the magnetic field from the 



center of the cloud by Ohmic dissipation. iNakano et al. I (120021 ) showed that the magnetic 



field is largely removed by Ohmic dissipation for 10 12 < n c < 10 16 cm -3 . Since the first core 
forms at n c ~ 10 12 cm -3 in this model, the first core is composed of the gas from which the 
magnetic field is already removed. Therefore, the first core has a large plasma beta (or weak 
magnetic field strength), as shown in Figure [T]c. 
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Figures[T](i and e show the density and velocity distributions around the first core at n c = 
6 x 10 20 cm -3 (second collapse phase). The cloud can collapse again for 10 16 cm" 3 < n r < 



10 c m 3 owing to dissociation of molecular hydrogen (i.e., second collapse; iMasunaga fc Inutsuka 



20001 ). The arrows in Figure [T]d and e indicate that the cloud collapses intensely around its 
center. The size of first core at this epoch is slightly smaller than that at the epoch of Fig- 
ure [TJa and b. However, the first core has a nearly spherical shape (e Q b — 1 at n c = 6 x 10 20 in 
Fig. [2]a) because the anisotropic Lorentz and centrifugal forces are weaker than the isotropic 
forces of the gravitational and thermal pressure gradient forces in this phase. The Lorentz 
force becomes weak due to magnetic dissipation, and an initially small rotation rate produces 
a weaker centrifugal force. 

Figure shows the angular velocity normalized by the free-fall timescale at the center 
of the cloud, whose val ue corresponds to the ratio of the rotational to gravitational energy 
( Machida et al. Il2005al ). When this value reaches Vt c / (AirGpc) 1 ^ 2 ~ 0.2, the rotational energy 



becomes c omparable to the gravitation a l energy, and cloud rotation b egins to affect cloud 



evolution (IMatsumoto fc Hanawal 120031 ; iMachida et al. I l2005al . l2006al ). In model SR, the 
effect of cloud rotation is significantly small in the second collapse phase (10 16 cm -3 < n c < 
10 20 cm -3 ) because Vt c / (AnGpc) 1 ! 2 <C 0.2. The magnetic field is again important for cloud 
evolution after the second collapse phase because Ohmic dissipation becomes less effective, 
and the magnetic field can be amplified. Figure [2]c shows that B c j ' p c has a minimum at 
n c ~ 10 15 cm~ 3 , and then continues to increase for n c > 10 15 cm~ 3 . However, the plasma 
beta inside the first core reaches only (3 V ~ 100 at n c ~ 10 21 cm -3 (FigJTJf), and therefore, 
the magnetic energy is much smaller than the thermal energy, even in the second collapse 
phase. 

Figures [T^ and h show the density and velocity distributions at n c = 2.6 x 10 21 cm" 3 
(the protostar formation epoch). In our calculations, the first core (i. e., shock layer) does not 



disap pear until the calculation ends. The second core (i.e., protostar; IMasunaga fc Inutsuka 



20001 ) forms inside the first core in Figures [Jj? and h. Since both the magnetic field and 
angular velocity continue to increase in the second collapse phase, the first core is slightly 
flattened both by Lorentz and centrifugal forces at the protostar formation epoch, as shown 
in Figure dpi. However, the protostar has an oblateness e Q b = 1.2 (Fig. [2]a) at its formation 
epoch (n c = 10 21 cm -3 ), and therefore, has a nearly spherical shape. The magnetic field 
and rotation period reach B = 2.18 kG and P = 3 days at the protostar formation epoch. 
The plasma beta around the protostar (or inside the first core) becomes /3 P ~ 10 — 10 3 
(Fig. []J). After the protostar formation (n c > 10 21 cm -3 ), the plasma beta becomes large 
(/3 P ~ 10 3 — 10 4 ) because the thermal energy becomes large [see eq. (J5J)]. The angular 
velocity normalized by the free-fall timescale, however, reaches fl c / 1 (AnGpc) 1 / 2 ~ 0.2 at the 
protostar formation epoch (Fig.[2]d); then, the cloud begins to rotate rapidly, and the rotation 
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significantly affects protostar evolution. 



3.1.2. Jet Formation 

Figure [3] shows cloud evolution after the protostar formation epoch (n c > 10 21 cm -3 ). 
Figures [3]a-c show the protostar (i.e., second core) surrounded by the shock layer as a dotted 
line. Figure [3]a shows the density and velocity distribution around the protostar 170 days 
after the protostar formation epoch. The protostar has a mass of 5.1 x 10 -4 M Q and a radius 
of R ~ 0.75 Rq at its formation epoch (n c ~ 10 21 cm -3 ). The arrows in Figure [3]a show 
that the gas around the protostar accretes spherically onto the protostar. Figure [3]& shows 
the protostar 207 days after its formation epoch. Figures [3]a and b show that the protostar 
becomes flat with time because the magnetic field and cloud rotation are amplified and affect 
protostar evolution, as shown in Figures [2]c and d. Two shocks (outer and inner) are seen 
in Figure [3]&. The outer shocks are located at z ~ ±6 x 10 _3 AU, while the inner shocks 
corresponding to the protostar are located at z ~ ±2 x 10~ 3 AU. As shown in Figure [3]c, the 
accretion speed inside the outer shock (\v\ ~ 20kms _1 in the region of z < |2 x 10 3 1 AU) is 
slower than that outside the outer shock (\v\ ~ 5 km s" 1 in the region of z > |2 x 10 3 1 AU). 
This structure of the shock and dist ribution of the de nsity and velocity, shown in Figure [3f), is 



similar to "the magnetic bubble" in iTomisakal (120021 ) . He showed that the magnetic pressure 



is amplified by the twisted magnetic field lines, and a bubble-like structure is formed in the 
weakly magnetized collapsing cloud. 

In model SR, the jet appears 215 days after the protostar formation epoch 0. Figure [3]c 
shows that the gas flows out from the central region along the vertical axis. The red lines in 
Figures W)-f indicate the border between the inflow and the outflow. Inside the red lines, the 
gas flows out from the central region [y z ^ for z ^ 0), and outside the redlines, it flows into 
the central region (v z ^ for z ^ 0) . The jet has a maximum speed of Vj et — 17 km s _1 at the 
same epoch as in Figure [3]c. The jet has an elongated structure and height-to-radius ratio 
of H/R « 3 in Figure [3]c. The height-to-radius ratio increases with time. Figures EH and 
e show the jet extending in the vertical direction keeping almost the same horizontal scale 
(R = 0.012 AU in Fig. [3]d and 0.018 AU in Fig. [3]e). The jet extends up to H = 0.043 AU 
with a maximum speed of fj e t — 22kms~ 1 at t c = 229 days (Fig. [3]d), while it extends up 
to H = 0.1 AU with Uj et ~ 26kms~ 1 at t c = 226 days (Fig. [3]e), where t c is the elapsed time 
from the protostar formation epoch (n c = 10 21 cm -3 ). Figure [3]f shows the structure of the 



2 For convenience, this paper refers to the flow driven from the first core as "the outflow" , and the flow 
driven from the second core (or protostar) as "the jet." 
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jet at the end of the calculation. Note the butterfly-like density distribution (white-density 
contour), which is caused by the strong mass ejection from the protostar. The jet has a well- 
collimated structure at this epoch. The height-to-radius ratio of the jet reaches H/R ~ 10. 
In Figure [3]f, the jet penetrates the first core, shown as a thick broken line, and extends up 
to the outside of the first core. At the end of the calculation, the jet reaches 0.2 AU with a 
maximum speed of t>j et ~ 32kms" 1 . 

Figure |H shows the shape of the jet (purple iso- velocity surface) and magnetic field lines 
(stream lines) at the same epoch as Figure[3]fbut with a different grid size (left: I = 22, right: 
I = 24). Figure H] left panel confirms that the jet has a well-collimated structure showing 
that the jet is strongly coiled by the magnetic field lines, with the bow shocks clearly visible 
near the upper and lower boundaries. Figure [3] (right panel) provides a close-up view of 
the left panel, indicating that the jet is driven from the protostar, which is shown by a red 
iso-density surface. We can see ripping den sity contours on the projected wall caused by 
strong mass ejection (IMachida et al. Il2006bl ). 

At the end of the calculation, the protostar in model SR has a mass of 4.04 x 10~ 3 M and 
a radius of 1.07 R Q , while the mass flowing out from the protostar is Mj et = 5.23 x 10~ 5 M & . 
At the same epoch, the protostar and jet have the angular momenta of J core = 4.23 x 
10 47 gcm 2 s _1 and Jj e t = 1-55 x 10 46 gcm 2 s _1 , respectively. Therefore, the specific angular 
momenta of the protostar and jet are j core = 5.25 x 10 16 cm 2 s _1 and jj et = 1.49 x 10 17 cm -3 , 
respectivel y This indicate s that the jet largely removes the angular momentum from the 
protostar. IMachida et al. I (120071 ) showed that, at its formation epoch, the protostar has a 
short rotation period compared with the observation. The jet may be one of the mechanisms 
for the angular momentum transfer after protostar formation. 



3.2. Cloud Evolution with Moderate Rotation 

3.2.1. First Core and Protostar Formation 

Figure [5] shows the cloud evolution for model MR after first core formation. Model 
MR has a parameter of uj = 0.03. The ratio of the rotational to gravitational energy for 
model MR is (3q = 3 x 10~ 3 , which is 100 times larger than that for model SR at its initial 
state (see Tabled]). In model MR, the first core forms at n c ~ 10 14 cm -3 , and has a mass 
of M ~ 0.01 Mq and a radius of R ~ 1 AU. Figures [5]a and [T]a show that the first core 
of model MR is larger than that of model SR, because the centrifugal force for model MR 
is stronger. The angular velocity, normalized by free-fall timescale at the center of the 
cloud, saturates fl c / (47rGp c ) 1/2 ~ 0.2 for n c > 10 14 cuT 3 (Fig.[2]d). This indicates that cloud 
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rotation significantly affects cloud evolution after the central density reaches n c ~ 10 cm" . 
Due to the rotation (or the centrifugal force), the first core for model SR, has a more oblate 
structure than that for model SR, as shown in Figure O?. The evolution of the oblateness for 
model MR is almost the same as that for model SR for n r < 10 14 cm~ 3 . On the other hand, 



the oblateness increases after the central density reaches n c 



10 15 cm 3 , 



while it maintains 



e b — 1 for n c > 10 14 cm -3 in model SR. Figure [5]c shows the plasma beta around the first 
core at n c = 3.2 x 10 16 cm -3 . The plasma beta for model MR is as large as that for model 
SR (/3p = 10 — 10 3 ) at the first core formation epoch (Fig. 0c). Thus, the magnetic field 
slightly affects cloud evolution inside the first core at this epoch. 



Figures [5]d and e show the density and velocity distribution at n c = 1.2 x 10 18 cm -3 
(the second collapse phase) . The arrows in Figure [BJd show that the cloud collapses while 
rotating: however, it collapses only in the radial direction for model SR (Fig. []]). The 
azimuthal component of the velocity (v^) inside the first core is comparable to or larger than 
the radial one (v r ) in model MR, while the azimuthal component of the velocity is much 
smaller than the radial one in model SR. Since the gas cloud inside the first core rotates 
rapidly, the central region increases its oblateness as the cloud collapses (Fig.[2]a). Therefore, 
the central region h as a m ore oblate or disk-like structure, similar to the first core shown in 



Saigo fc Tomisakal (120061 ). in which they investigated the equilibrium state of the rotating 



first cores. In general, since the angular momentu m is transferred by magnetic braking in 
the magnetized clouds (IBasu fc Mouschoviaslll994l ). the cloud rotation slows with time. In 
this model, however, the cloud rotation hardly slows because the magnetic energy around 
the first core is too small to effectively work the magnetic braking. The small energy of 
the magnetic field is realized by Ohmic dissipation. Therefore, the cloud around the first 
core has a large plasma beta (/? p > 10; Fig. [5]c). Figure [2]c compares the magnetic fields 

1 /2 

normalized by the density at the center of the cloud (B c / p c ) for model MR (broken line), 
which are smaller than those for model SR (solid line) because for model MR, the cloud 
collapses, maintaining the disk-like structure, while that in model SR collapses spherically. 
The growth rate of the magnetic field depends on the geom etry of the collapse and is muc h 



larger for a spherical collapse than for a disk-like collapse (iMachida et al. Il2005al . l2006al ). 
As a result, the protostar has a magnetic field of B = 0.4 kG for model MR which is five 
times smaller than that for model SR at its formation epoch (n c 



HP cm" 3 ). 

Figures \5\g and h show the density and velocity distribution at n c = 2.3 x 10 21 cm -3 
(the protostar formation epoch). In each panel, the protostar (i.e., the second core) is rep- 
resented by a red-dotted line. In model MR, the protostar forms at n c ~ 2 x 10 20 cm~ 3 , 
and has a mass of M ~ 3.6 x 10~ 4 M Q and a radius of R = 1.5 R & . As shown in Fig- 
ure [5fi, at the protostar formation epoc h, the first core sags dow nward in the center in a 
concave structure similar to model C in ISaigo fc Tomisakal (120061 ). This concave structure 
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is considered to be caused by rapid rota tion and an increase in the density of the first core 
(for details, see lSaigo fc Tomisakall2006l ). As shown in Figure [5ft, since the gas density falls 
below n c < 10 12 cm -3 in the regions just above and below the first core at the center of the 
cloud (\x\ < 0.1 AU), Ohmic dissipation is less effective i n these regions. Note that Ohmic 
dissipation is effective for 10 12 cm -3 < n c < 10 15 cm" " 3 jNakano et al. Ibooj ). Therefore, 
the magnetic field is strong around the protostar (/3 P ~ 0.1). Since the central region be- 
comes increasingly thinner and has an increasingly concave structure over time, the magnetic 
field around the protostar barely dissipates through Ohmic dissipation. Therefore in model 
MR, the torsional alfven wave generated by the rotation of the protostar can be transferred 
outside without dissipation. In model SR, however, the magnetic field dissipates and the 
plasma beta maintains a high value (/3 P ~ 10 — 10 3 ) around the protostar, because the cloud 
collapses spherically with slow rotation, and Ohmic dissipation is effective around the proto- 
star. Figure [5ft shows that the outflows are driven near the first core. The region where the 
outflow appears in Figure [5ft corresponds to the low-plasma beta region (/3 P < 0.1) shown in 
Figure O. 



3.2.2. Outflow and Jet formation 



Figure M shows the evolution of outflow for model MR. Figures Eft-c show the density 
and velocity distribution around the first core on the y = plane. The gas accretes onto the 
first core (white dotted lines) in the whole region outside the first core (Fig.[6]a), while the gas 
in the regions above and below the first core (inside the red lines) outflows from the central 
region (Fig. [6ft) . In the ideal MHD regime, the outflow in the collapsing cloud appears in close 



proximity to the first core ([To misaka 120021 ; iMatsumoto fc Tomisakal 12004 ; iMachida et al. 



2004j . l2005bl ; lBanerjee fc Pudrital2006l ). However, as shown in Figure Eft, the outflow appears 
far from the first core in the non-ideal MHD regime because the magnetic field near the first 
core is not well coupled with the neutral gas due to Ohmic dissipation. Thus, the torsional 
alfven wave caused by the rotation of the first core is not effectively transferred. However, 
the torsional alfven wave generated far from the first core can maintain itself, because the gas 
density is low and Ohmic dissipation is less effective in these regions. As the cloud collapses, 
the outflow can be driven even near the first core, as shown in Figure [6]c, because in these 
regions, the gas density gradually decreases and Ohmic dissipation becomes less effective 
over time, as shown in Figures Eft, e, and h. Finally, the outflow is anchored by the first 
core, as seen in the ideal MHD calculations. At the end of the calculation, the outflow has a 
maximum speed of 3.2 kms" 1 , which is comparable to the results of ideal MHD calculations 
(lTomisakall2002l ) . Therefore, the outflow in our calculation (the non-ideal MHD calculation) 
differs from that of the ideal MHD calculation in the early phase, while the outflow features 



-14- 



are almost the same (e.g., the shape and speed of the outflow) for both the non-ideal and 
ideal MHD calculations in the later phase. 

Figure [7J shows the evolution of a jet driven from the protostar. In this figure, the first 
core and protostar (i.e., the second core) are shown as dotted white and red lines, respectively. 
The protostar has a disk-like structure, while the first core has a concave shape (Fig. [71a). 
The gas flows out from the regions above and below the protostar (Fig. [7J&). Figure [73 shows 
horn-like structures at z = ±0.1 AU, which are caused by strong mass ejection. This panel 
shows that only the gas located in the regions just above and below the protostars flows out 
from the protostar, while a large fraction of the gas accretes onto the protostar. At the end 
of the calculation, a strong jet with a maximum velocity of 17.2 km s" 1 is driven from the 
protostar, as shown in Figure [TJc. 

Figures [S]a-c (upper panels) show the density (color and contours) and velocity distri- 
butions (arrows) around the protostar for model MR. The velocity of ^-component (v z ; color 
and contours) and magnetic field (arrows) corresponding to each upper panel of Figure [8] are 
plotted in Figures E^/ (lower panels). The first core and protostar (i.e., second core) are 
plotted as dotted white and red lines, respectively. The outflow extends from the region just 
above the first core to the outside of the I = 17 grid boundary (Fig. [8]a). In model MR, the 
outflow extends to 7 AU at the end of the calculation. 

There are two velocity peaks in Figure EK The upper one (peak 1) has a velocity 
of v z ~ 2.5 kms -1 at (x,z) ~ (±0.7 AU, 5AU), while the lower (peak 2) has a velocity 
of v z ~ 3kms _1 at (x, z) ~ (±0.2 AU, 1AU). Figures [8]6 and e provide close-up views of 
Figures [8]a and d. In Figure [8]e, strong flows (peak 3) appear near the surface of the first core 
[(x,z) ~ (±0.1 AU, 0.3 AU)]. Figures [8]c and / show close-up views of Figures [8f> and e. A 
strong jet is driven from the region above the protostar in Figure [H]c, corresponding to peak 
4 in Figure [S]f. This peak (peak 4) has a velocity of v z ~ lSkms" 1 at (x, z) ~ (±0.02 AU, 
0.08 AU). The lower panels of Figure [8] show four velocity peaks (peak 1-peak 4). These 
peaks imply that speeds (strength) of the outflow and jet change with time, because they 
represent a history of the outflow and jet. We confirmed that peaks 1 and 2 appear before 
protostar formation, and peaks 3 and 4 appear after protostar formation. We also confirmed 
that each peak is related to the oscillation of each core. For example, peaks 1 and 2 appear 
every time the first core oscillates. It is, therefore, considered that the outer peaks (peaks 
1 and 2) originate from the first core, while the inner peaks (peaks 3 and 4) originate from 
the protostar. Thus, in our definition, the outer peaks correspond to the outflow, while the 
inner peaks correspond to the jet. 

Figure [9] shows the structure of the outflow and jet, and the configuration of the magnetic 
field lines. It also shows the shapes of the first core (left panel; the projected density 
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contours on the wall) and the protostar (right panel; the red isosurface). The purple and 
blue surfaces in Figure [9] indicate the iso- velocity surface of v z — 5 km s _1 and v z = 0.5 km s _1 , 
respectively. The flow inside the purple iso-velocity surface has a velocity of v > 0.5 km s" 1 
(low- velocity component; LVC), while the flow inside the blue iso-velocity surface has a 
velocity of v > 5 km s" 1 (high- velocity component; HVC). The HVC is enclosed by the LVC. 
The LVC flow is mainly driven from the first core, while the HVC flow is mainly driven from 
the protostar. The LVC and HVC are strongly coiled by the magnetic field lines anchored 
to the first core and protostar, respectively. 

In model MR, the protostar has a mass of 4.52 xlO -3 M and a radius of 1.07 _R , while at 
the end of the calculation, the mass flowing out from the protostar is M out = 1.42 x 10~ 3 M & . 
At the same epoch, the protostar and outflow gas (outflow and jet) have angular momenta 
of J C ore = 1-53 x 10 48 gcm 2 s _1 and J out = 1.93 x 10 48 gcm 2 s _1 , respectively. Therefore, the 
specific angular momenta of the protostar and outflow gas (i.e., the outflow and jet) are 
jcore = 1-67 x 10 17 cm 2 s _1 and j ont = 6.82 x 10 17 cm -3 , respectively. 

3.3. Cloud Evolution with Rapid Rotation 

3.3.1. First Core Formation 



Figure [10] shows the cloud evolution for model FR after the first core formation. Model 
FR has a parameter of uj = 0.3. The ratio of the rotational to gravitational energy for 
model FR is (3q = 0.3, which is 10 4 times larger than that for the initial stage of model 
SR(Table[T|). The first core forms at n c ~ 8 x 10 12 cm -3 , and has a mass of M ~ 0.026 M & 
and a radius of R ~ 5.68 AU at its formation epoch. The mass and size of the first core 
for model FR are larger than those for models SR and MR because the cloud is rotating 
rapidly. The cloud in model FR has VL C / (AnGpc) 1 ! 2 ~ 0.2 in the initial stage, and maintains 
this value for n c < 10 12 cm -3 , as shown in Figure [2]g?, which indicates that the centrifugal 
force significantly affects the formation of the first core. Owing to rapid rotation, the cloud 
has a disk-like structure near the center, as shown in Figure [Tol a. The oblateness of model 
FR increases rapidly, as shown by the dotted line in Figure [2]a, and saturates e ^ — 5 — 6 
for n c ~ 10 7 cm -3 . In this model, the first core has a disk-like shape at its formation epoch. 
Then, the disk-like first core deforms to a ring, an d the ring fragments into several pieces at 
n c ~ 3 x 10 14 cm -3 . iMachida et al. I (120041 . l2005bl ) showed that fragmentation occurs in the 
rapidly rotating cloud. Since we focus on the mechanism of the outflow and jet, we stopped 
the calculation when fragmentation occurred. We will discuss the fragmentation process of 
the first and second cores in a subsequent paper. 
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3.3.2. Outflow formation 

Figure [TOfe shows the density and velocity distribution around the center of the cloud 
330.3 yr after the first core formation. This panel shows that outflow is driven from the center 
of the cloud. The outflow extends up to 120 AU with a maximum speed of i> ou t — 2.5 km s -1 
at this epoch. Figure fTUb shows the structure of the outflow 476.1 yr after the first core 
formation. At the end of the calculation, the outflow extends up to 200 AU with a maximum 
speed of v out = 3.1 kms" 1 . Figures [TDl). [3j and [8] show that the opening angle of the outflow 
for model FR is larger than that of the jet for models SR and MR. The height-to-radius 
ratio of the outflow is H/R ~ 2.2 (Fig. UM>) to 2.5 (Fig. [He) in model FR, while that 
of the jet is H/R ~ 10 in model SR. In model FR, the outflow maintains its shape (i.e., 
H/R ~ constant) regardless of time, while the jet in model SR becomes slender along the 
vertical axis with time (i.e., H/R increases with time). 

Figure [H] shows the density and velocity distribution on the y = (upper panel) and 
z = (lower panel) planes at the same epoch as Figure [TUb . but with different grid scales. 
Figure [Tib is 16 times magnification of Figure [TTr and shows the ripping density contours 
caused by the strong mass ejection from the center of the cloud. Since the outflow interrupts 
gas accretion in the regions above and below the central core, the gas accretes onto the 
central region only via a thin disk (Fig. ITTb). The upper panel in Figure [TTb shows that 
the outflow is driven from the first core (dotted line). At this epoch, the first core, which 
has a disk-like shape at the formation epoch, has deformed to a ring-like shape (Fig. HTb 
lower panel). As shown in Figure [Tib, the density of the first core exceeds n c > 10 12 cm -3 , 
while the density of the ambient gas outside the first core is n r < 10 12 cm ~ 3 . Since Ohmic 



dissipation is effective for 10 12 cm -3 < n c < 10 15 cm -3 (INakano et al. 1 120021 ) . the dissipation 
of the magnetic field becomes effective inside the first core, while it is barely dissipated by 
Ohmic dissipation outside the first core. Therefore, the magnetic field lines can be twisted 
just outside the first core, and the torsional Alfven wave transfers to the outside, as shown 



in the ideal MHD calculations (ITomisakal 120021 ; iBanerjee fc Pudritzll2006l ). 



Figure [T21 shows the shape of the outflow (blue iso- velocity surface) and the configuration 
of the magnetic field lines (streamlines) at the same epoch as Figure [TTJ This figure shows 
that the magnetic field lines are strongly twisted inside the outflow, while they are loosely 
twisted outside the outflow. Due to rapid rotation, the first core forms at a lower density 
for model FR (n c ~ 8 x 10 12 cm -3 ) than for models SR and MR (n c ~ 10 14 cm -3 ); therefore, 
Ohmic dissipation is less effective outside the first core in model FR. The magnetic field lines 
anchored to the rapidly rotating first core are effectively twisted, and the Lorentz force can 
drive the strong outflow near the first core. We stopped the calculation for model FR at 
n c ~ 10 15 cm -3 , because fragmentation occurred. Thus, we cannot determine further cloud 



-17- 



evolution for model FR. We expect that each fragment has a small spin angular momentum 
because of redistri bution of the angular momentum (spin and orbital angular momenta) 
for fragmentation (jMachida et al. 2005b ). and the magnetic field is dissipated by Ohmic 



dissipation in each fragment. Thus, we expect that each fragment traces a similar evolution 
to models SR and MR. 

In model FR, the first core has a mass of 9.42 x 10~ 2 M and a radius of 5.32 AU, 
while the mass flowing out from the first core is M out = 6.01 x 10~ 3 M at the end of the 
calculation. At the same epoch, the first core and outflow have the angular momenta of 
■/core = 8.48 x 10 51 gcm 2 s _1 and J out = 4.95 x 10 50 gcm 2 s~ 1 , respectively. Therefore, the 
specific angular momenta of the first core and outflow are j com = 4.52 x 10 19 cm 2 s _1 , and 
Jout = 4.14 x 10 19 cm -3 , respectively. Thus, in model FR, the outflow largely removes the 
angular momentum from the first core. 



3.4. Relationship between Outflow / Jet and Initial Cloud Rotation 



In the previous sections, we showed the evolutions of three di fferent clouds, while we 

2005al . 



investigated the cloud evolutions of 80 models in total. As shown in Machida et al. { 



2006al . 120071 ). cloud evolution depends only on the initial angular velocity in the magnetic- 
force dominant clouds Vt Q / B Q < O cri /5 cri = 0.39 G 1 / 2 cj 1 , where VL Q , B , and c s are the initial 
angular velocity, magnetic flux density, and sound spee d, respectively. Our choice o f these 
strongly magnetized clouds is supported by observations (jCrutcherlll999l ; ICasellill2002l ) . From 
our calculations, we found the emergence conditions for the outflow and the jet fall into three 
categories, depending on the initial rotation rate uj: 



(1) uj > 0.01, both outflow and jet can be driven from each core, 

(2) 0.002 < uj < 0.01, only a jet can be driven from the protostar, and 

(3) uj < 0.002, neither outflow nor jet can be driven from each core. 



The outflow always appears when condition (1) is satisfied. However, the jet does not 
appear in some models even when condition (1) or (2) is satisfied. Thus, condition (1) is 
the necessary and sufficient condition for driving outflow, while condition (1) and (2) is the 
necessary condition for driving jet. For example, although 24 models satisfy condition (2), 
the jet appeared only in nine models. Although we followed the evolution of the accreting 
protostar up to ~200 CPU hours for each model, we did not observe any sign of the jet 
in those nine models. After the protostar is formed, the magnetic field becomes strong 
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in the regions above the poles of the protostar; thus, the plasma beta becomes ^ < 1 
in these regions. Therefore, it is difficult to calculate for a long period after protostar 
formation because the time step becomes significantly short owing to the increased Alfven 
speed. Therefore, we could not determine the reason the jet does not appear in some models. 
We could not determine if the jet would never appear in those models or if it might appear 
much later. To clarify the long-term evolutions, we need to develop a calculation method 
with implicit time integration. 



4. Discussion 

4.1. Driving Mechanism of Outflow and Jet 

In this paper, for convenience, we called the flows driven from the first core the outflow, 
and the flows driven from the protostar (i.e., the second core) the jet. In model SR, no outflow 
appears because the initial cloud rotates slowly, and the rotation rate is not sufficiently 
amplified at the first core formation epoch. On the other hand, we cannot observe the 
jet driven from the protostar in model FR because we stopped calculation once it showed 
fragmentation. In model MR, both the outflow and the jet appear. In this section, to 
investigate the driving mechanism of the outflow and the jet, we focus on the flows that 
appeared in models SR (the jet) and FR (the outflow). We do not discuss the flows in model 
MR because it is difficult to separate the jet from the outflow as shown in Figures EHE 

At the end of the calculation, the jet has a well-collimated structure with high speed 
(vj e t ~ 30kms _1 ) and a large height-to-radius ratio of H/r > 20 . The H/r ratio of the jet 
increases with time, as shown in §3.1.2 and 3.2.2. The outflow, however, has a wide opening 
angle with slow speed (t> out — 3kms _1 ), and expands outwardly, maintaining the height-to- 
radius ratio of H/r ~ 2 — 2.5. The difference of the speeds between the jet and outflow can be 
understood from the difference in the gravitational potential. In model SR, the protostar has 
a mass of M ~ 4.04 x 10~ 3 M & and is a radius of R ~ 1.07 i? at the end of the calculation. 
The Kepler speed corresponding to the mass and radius is fRepier = 26kms _1 , which is 
comparable to the speed of the jet in model SR (vj et — 30kms _1 ). On the other hand, at 
the end of the calculation for model FR, the first core has a mass of M ~ 9.42 x 10~ 2 M 
and a radius of R ~ 5.32 AU. Thus, the Kepler speed of the first core is ^Kepler = 3.96 km s _1 , 
which is comparable to the speed of the outflow (v ont — 3kms _1 ). As a result, both the jet 
and the outflow have flow speeds similar to the Kepler speed of their respective cores. 
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4-1.1. Jet Driven from the Protostar 

The difference in the degree of collimation between the outflow and the jet (the well- 
collimated structure of the jet and the wide opening angle of the outflow) can be understood 
from the difference in their driving mechanisms and the configuration of the magnetic field 
lines around their drivers (i.e., the fist core and protostar). Figure [T31 shows the shapes of the 
jet (model SR; left panels) and outflow (model FR; outflows) with thick red lines indicating 
the border between the accretion and the jet or outflow. This shows that the collimation of 
the jet (left panel) is stronger than that of the outflow (right panel). 

First, we discuss the collimation and driving mechanism of the jet. The upper panel 
of Figure [13] show the ratio of the toroidal to poloidal components of the magnetic field 
(£?toroidai/-B P oioidai) at each mesh point around the outflow (model SR; left) and the jet 
(model MR; right). In Figure [TBla. the toroidal field dominates the poloidal field inside the jet 
(^toroidal/ ^oioidai - 10), while the toroidal field barely exists outside the jet (-Steroidal/ Ppdoidai ~ 
0.01). In model SR, the magnetic field is dissipated by Ohmic dissipation, and decoupled 
from the neutral gas inside the first core (or outside the second core), as shown in Figure [ljj. 
Thus, the magnetic field lines are relaxed by the magnetic tension force and are almost 
straight [B z 3> B r , in the cylindrical coordinates), because magnetic field lines move 
freely, irrespective of the neutral gas. Even after the magnetic field is well-coupled with 
the neutral gas for n c > 10 15 cm -3 , the straight magnetic field lines are distributed around 
the protostar, which reflects past Ohmic dissipation. Therefore, the vertical component of 
the magnetic field (B z ) is dominant outside the protostar. The lower panels in Figure [TBI 
show the plasma beta around the jet (left; model MR) and the outflow (right; model FR). 
Figure [T3J6 shows that the magnetic field is weak (/3 P ~ 10) outside the jet (outside the red 
line), while it is very strong (j3 p < 0.01) inside the jet (inside the red line), compared to the 
thermal pressure. The sharing motion between the protostar and ambient gas amplifies the 
magnetic field inside the jet, since the formed protostar rapidly rotates because the mag- 
netic braking is less effective, explained in §3.1. In general, when the magnetic field is weak 
around the driver, it is considered that the magnetic pr essure gradient force dominates the 



magnetocentrifugal force in terms of driving the flow (jUchida Sz Shibata I Il985l ; iTomisaka 



2002( 1. 



To investigate the driving mechanism of the jet and the outflow, we calculate the Lorentz 
(F Lor entz), centrifugal (-^centrifugal), and thermal pressure gradient (-^pressure) forces parallel to 
the poloidal component of the magnetic field (-Bp). When the self-gravity term in equation (jSj) 
is ignored, the equation of motion can be rewritten as 

p-^ = -p(v ■ V) - VP - - [B x („ x B)} . (7) 
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We consider the dot product between each term in the left-hand side in the equation (j7]) and 
the poloidal component of the magnetic field (B p ). Thus, each force parallel to the B p is 
written as 

^centrifugal = \p(v • V) • e p | , (8) 
-^thermal = |VP-e p |, (9) 

F Lo rentz = — \B X (v X B) ■ e p \, (10) 
47T 

where e p = B p /\B p \. In the left panels of Figure [TH we plot the ratio of the Lorentz 
to the centrifugal force (FLorentz /F cen trifugaij Fig- EK), and the ratio of the Lorentz to the 
thermal pressure gradient force (-^Lorentz/ -^thermal; Fig. [life ) in the outflow region at each 
mesh point, where the outflow region is defined as meshes with v z > O.lkms -1 for z > 0. 
The right panels (Fig. [TJ]c and d) show the plots for the same ratios of F Lorentz / F ccntrifngSb i, 
and -FLorentz /thermal in the accreting region (v z < 0.1 km s^ 1 for z > 0). In Figure [T4"l red 
dots and lines are data from the region around the jet in model SR, while black dots and 
lines are data from the region around the outflow in model FR. In Figure dU left panels 
show that, inside the jet (red dots and line), the Lorentz force is equivalent to the centrifugal 
force ( FLorentz /-^centrifugal — 1; Fig. [Tib), while the Lorentz force is about ~ 30 — 40 times 
stronger than the thermal pressure gradient force (-FLorcntz/Fccntrifugai — 30 — 40; Fig. [life) . In 
the accreting region around the jet (red dots and lines in Fig. [T4b and d), the Lorentz force 
is much weaker than both the centrifugal (-FLorentz /-^centrifugal — 10~ 5 — 10~ 4 ) and thermal 
pressure gradient forces (FLorentz /F P ressure — 10~ 6 — 10~ 4 ). Thus, the Lorentz force is not 
effective in the accreting regions. 

In model SR, a rapidly rotating protostar forms because the magnetic field around the 
center of the cloud is weak and magnetic braking is less effective. Thus, the magnetic field 
lines anchored to the protostar make a strong toroidal field for the rapid rotation of the 
protostar. The toroidal field caused by the protostar transmits in the vertical direction as 
the torsional alfven wave, generating a large magnetic pressure gradient along the vertical 
axis. The jet is considered to be driven by this strong magnetic pressure gradient force. In 
general, the jet has good collimation when the magnetic pressure gradient force is domi- 
nant (-Btoroidai 3> -Bpoioidai) for drivi ng the flow , beca use the magnetic field lines are pinched 



by the toroidal field, as shown by iTomisakal (120021 ) . Furthermore, the jet in model SR is 



guided by the straight configuration of the magnetic field lines (B z ^> B r , B p ) outside the 
jet region. Thus, the well-collimated jet is caused by both the driving mechanism of the 
magnetic pressure gradient and the straight configuration of the magnetic field lines around 
the protostar. 
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4-1.2. Outflow Driven from the First Core 



Figure [13] shows the ratio of the toroidal to poloidal components of the magnetic field 
( .Btoroidai /-Bpoioidai; Fig- HSfc) and the plasma beta (/3 P ; Fig. fTBIcf) around the outflow for model 
FR. The toroidal field is ~ 3 — 5 times stronger than the poloidal field (-Btoroidai/-B po ioidai — 
3 — 5) inside the outflow (inside the red line). However, this ratio is smaller than that of 
the jet ( -B toro idai /-Bpoioidai — 10; model SR) because the rotation rate in model FR, which 
causes the toroidal field, is small around the first core, due to effective magnetic braking. 
Outside the outflow, however, the poloidal field is slightly larger than the toroidal field 
( -Steroidal /-Bpoioidai — 0.3). Figure [T3lc shows that the magnetic energy is larger than the 
thermal energy inside the outflow (/3 P < 0.01), while the magnetic energy is comparable 



to or slightly weaker than the thermal energy outside the outflow (f3 v 



0.1 - 1). The 



gas around the first core has a density range of 10 8 cm 3 < n c < 10 12 cm 3 in model SR 



(Fig. [TO]) . Since Ohmic dissipation is ineffective in thi s region, the magnetic energy b ecomes 
comparable to the thermal energy as demonstrated by iMachida et al. I (j2005al . l2006al ) . After 
the first core formation, the toroidal field can be amplified because the rotational timescale is 
shorter t han the collapsing time scale inside the first core, and the magnetic field lines begin 
to twist (IMachida et al. Il2006dl ). Therefore, the magnetic energy inside the outflow (or near 
the first core) becomes larger than that outside the outflow (or far from the first core). 



thermal ^ 



As shown in Figures [TH a and b, the centrifugal force is dominant over the Lorentz and 
thermal pressure gradient forces inside the outflow (-^Lorentz /-^centrifugal < 1, and Fi, OTersk7 ,/ F\ 
1; black dots and lines), while it is comparable to the Lorentz force inside the jet (red 
dots and lines). Outside the outflow, both the centrifugal and thermal pressure gradi- 
ent forces are stronger than the Lorentz force (i^entrifugai — -^thermal > ^Lorentz)) since 
^Lorentz /Centrifugal ^ 0.1 (Fig. [He) and F Lorentz /F th ermai ^ 0.1 (Fig. EH). These results 
indica te that the outflow in mo del SR is mainly driven by the magnetocentrifugal mecha- 
nism ( Blandford fc Payne! 1982 ). which causes the loosely pinched magnetic field lines inside 
the outflow. In addition, before the outflow appears, the magnetic field lines around the first 
core expand sideways toward the rotation axis from the c enter of the cloud, whose shape 
looks like a capital letter U or V, as shown by iTomisakal (120021 ) . Thus, the collimation of 
th e magnetic field l ines be comes worse with distance from the center, as shown in Figure 8a 
of IMachida et al. I (l2006dl ). Since the outflow traces this configuration of the magnetic field 
lines, it does not have good collimation and the outflow driven from the first core has a wide 
opening angle because of the magnetocentrifugal driving mechanism and the configuration 
of the magnetic field lines. The former realizes the loosely pinched magnetic field lines, and 
the latter swells the outflow. 
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4.2. Dependence on the Initial Magnetic field and Rotation Rate 



In this paper, we investigated cloud evolution with same magnetic field strengths (B^ = 
17 nG) but different rotation rates [fi ini = 7 x 10 -16 s -1 (model SR), 7 x lO^s -1 (model 
MR), and 7 x 10 _14 s _1 (model FR)]. As shown in previous section, the speeds of the jets 
and outflows depend on the mass and size of their drivers. The mass and size of the first 
and second cores at their formation epochs depend on the initial rotation rate. Therefore, 
the speeds of the jet and outflow depend on the initial rotation rate of the cloud core. 
For example, in a rapidly rotating clouds, the first core forms at lower density with larger 
size. When the first core has a large size and low density, the outflow is slow because the 
gravitational potential is shallow. On the other hand, a slowly rotating cloud produces a 
first core with smaller size and higher density. In this case, there is high-speed outflow, 
owing to the deep gravitational potential of the first core. Thus, the outflow has higher 
speed in a more slowly (rapidly) rotating cloud. In addition, a cloud with an extremely low 
rotation rate produces no outflow because the first core has such a small angular momentum 
at its formation epoch that the magnetic field lines are barely twisted. The speed of the jet, 
however, depends on the initial rotation rate only slightly because the angular momentum 
around the center of t he cloud converges to a certain value by the time the protostar (i.e., 
second core) is formed (IMachida et al. 1 120071 ). Thus, the variation caused by the rotation of 
initial cloud is not a factor affecting the time of the protostar formation epoch. Note that 
the properties of the jet and outflow might possibly change much later because we followed 
the evolution of these flows for only a short period. 

Does the difference of magnetic field st r ength i n the i nitial cloud affect the proper- 
ties of the jet and outflow? IMachida et al. I (l2005al . l2006al ) showed that the evolution of 
the molecular cloud depends only on the ratio of the angular velocity to the magnetic flux 
density of the initial cloud (f^m/Ani)- Initial differences in the distribution of the density, 
velocity, and magnetic field hardly affect cloud evolution, especially once the central region 
becomes adiabatic (n c > 10 11 cm -3 ). Cloud evolution is controlled by the magnetic field when 
foni/flini < 0.39 G 1/2 c7\ while it is controlled by the rotation when tt ini /B ini > 0.39 G l/2 c~ l 



(IMachida et al. 112006a ). Although we can choose arbitrary parameter of Vt m \ and B^, ob- 



servations indicat e that the magnetic energy i s much larger than the rotation energy in many 



molecular clouds (jCrutcherlll999l ; ICasellill2002T) . and clou ds haye ^ ini /i?i n i > 0.39 G 1 ^^ 1 (i.e 



the magnetic- force dominant clouds in IMachida et al. Il2005al ). In magnetic-force dominant 
clouds, the magnetic fiel d converges to a specific value in th e isothermal phase (the mag- 
netic flux - spin relation; IMachida et al. Il2005al . l2006al . 120071 ). Thus, cloud evolution does 
not depend much on the initial magnetic field strength. As a result, outflows and jets with 
similar properties may appear in these clouds when they have the same angular velocities at 
their initial state. 
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4.3. Duration of Driving Outflow 

In this paper, we have shown that the outflow is driven from the first core. After the first 
core formation epoch, the collapse timescale becomes longer than the rotational timescale, 
the magnetic field lines are strongly twisted for the rotating first core, and outflow appears 
as shown in §3.2 .2. Observations have shown that outflows have lifetimes of ~ 6 x 10 4 yr 



( TWu et al. II2004I ). However, since we followed cloud evolution only through the very early 
phase, we could not determine whether outflow continues to be driven for a long time. The 
spherical symmetric calculation indicates that the first core disappears ~ 1000 yr after its 
formation epoch. In model FR, we followed cloud evolution for ~ 800 yr after the first core 
is formed. In this model, the first core remains until the end of the calculation. The outflow 
that appears just after the first core formation epoch continues to be driven from the first 
core. The rotating first cores seem to have lifetimes of several thousand years or more (e.g., 



Saigo &: Tomisakal 120061 ) . Thus, outflows driven from these cores may maintain themselves 
for several thousand years. Note that it is expected that a jet will continue to be driven 
from the protostar once it appears because the protostar does not disappear. 

After the first core disappeared, we could not determine whether the outflow continues 



to be driven. ISaigo fc Tomisakal (120061 ) suggested that after the first core disappears, the 
remnant of the first core forms a disk-like structure or torus. Outflow may continue to be 
driven from these objects even after the disappearance of the first core. However, if the 
outflow does not continue after the first core disappears, then outflow, as shown in this 
paper, might be a transient phenomena, in which the outflow is driven only for ~ 1000 yr. 
In this case, the flow from the first core, which we call outflow in this paper, might be 
a special phenomenon that would be observed only in limited numbers of young stellar 
objects. To understand the duration of the outflow driven from the first core, further long- 
term calculations, possibly with an implicit code, are necessary. At present we conjecture 
that the outflow in our calculations will be driven for > 10 4 years. 



5. Summary 

Observation shows that the outflows have wide opening angles and low flow speeds, 
while jets have good collimation and high flow speeds. Outflow has been considered to be 
entrained by the jet driven from a circumstellar disk around the protostar. In this paper, 
we calculated the cloud evolution from the molecular cloud core to protostar formation. As 
a result of calculations, we found that two distinct flows (outflow and jet) are driven from 
different objects (the first, and second cores), and the features of outflow and jet (flow speed 
and collimation factor) were naturally reproduced. 
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Our results show that the flow appearing around the first core has a wide opening angle 
and slow speed, while the jet appearing around the protostar has a well-collimated structure 
and high speed, as shown in Figure [Tj3 The speed difference is caused by the difference 
of the depth in the gravitational potential. The flow speed (i.e., the speed of the outflow 
and jet) corresponds to the Kepler speed of each object (the first and the second cores). 
Because the first core has a shallow gravitational potential, its flow (outflow) is slower. 
The flow (jet) driven from the protostar, which has a deeper gravitational potential, has 
a high speed. In our calculations, the outflow and the jet have speeds of v out ~ 3kms _1 
and Vjet — 30kms _1 , respectively. These speeds are slower than those of observations. 
Typically, observed molecular outflow and optical jet have speeds of f ou t,obs — 30kms _1 , 
and fjet.obs — lOOkms -1 , respectively. However, since the first and second cores (protostar) 
have mass of M firstcore = 0.01 M Q and M secondcore ~ 1O _3 M , respectively, at the end of the 
calculations, each core increases its mass in the subsequent gas accretion phase. The Kepler 
speed increases with the square root of the mass of the central object at a fixed radius. When 
the mass of each core increases by 100 times, the Kepler speed increases 10 times. Thus, 
the speed of the outflow and jet may increase by 10 times, and reach t> out ~ 30kms _1 and 
ujet = 300 km s -1 , respectively, which correspond to typical observed values. 

The outflow has a wide opening angle, while the jet has a well-collimated structure. 
This is caused both by the configuration of the magnetic field lines around the drivers and 
their driving mechanisms. The magnetic field lines around the first core have an hourglass 
configuration because they converge to the cloud center as the cloud collapses, and Ohmic 
dissipation is ineffective before the first core formation. The flow appearing near the first core 
is mainly driven by the magnetocentrifugal wind mechanism (disk wind). The centrifugal 
force is dominant in the outflow, whereas near the protostar, the magnetic field lines are 
a straight, and the magnetic pressure gradient mechanism is more effective for driving the 
jet. The magnetic field lines are stretched by the magnetic tension force near the protostar 
because the magnetic field is decoupled from the neutral gas. However, the magnetic field 
lines are strongly twisted in the region in close proximity to the protostar, where the magnetic 
field is coupled with the neutral gas again. Thus, the strong toroidal field generated around 
the protostar can drive the jet, which is guided by the straight configuration of the magnetic 
field. 

Our calculations do not completely reject the well-known concept that the observed 
molecular outflow is entrained by the jet, because we calculate the formation of the jet and 
outflow only in the early star-formation phase. Further long-term calculations are needed to 
understand the mechanism of the jet and outflow in more detail. 

We have greatly benefited from the discussion with T. Nakano and K. Saigo. We also 
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Table 1: Model parameters and calculation results 



Model 




B [mG] 


Q [s- 1 ] 




M O [M ] 


a Q 


A) 


7o 


Bf (kG) 


P (day) 


SR 


0.003 


17 


7.0 x 10" 


-16 


7.6 


0.5 


3 x 10" 5 


0.96 


2.18 


3.0 


MR 


0.03 


17 


7.0 x 10" 


-15 


7.6 


0.5 


3 x 10" 3 


0.96 


0.40 


2.1 


FR 


0.3 


17 


7.0 x 10- 


-14 


7.6 


0.5 


3 x IQ" 1 


0.96 
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Fig. 1. — Time sequence of model SR. Upper panels: Density (color-scale) and velocity 
distribution (arrows) on the cross-section in the z = plane (a, d, and g). Middle panels: 
Density (color-scale) and velocity distribution (arrows) on the cross-section in the y = 
plane (b, e, and h). Lower panels: Plasma beta (color-scale and contours) and magnetic field 
(arrows) on the cross-section in the y = plane (c, /, and i). Panels from left to right are 
snapshots at the stages of n c = 6.5 x 10 16 cm -3 (I = 20, panels a, b, and c), 6.0 x 10 20 cm -3 
(I = 20; panels d, e, and f), and 2.6 x 10 21 cm~ 3 (/ = 20; panels g, h, and i), respectively, 
where I denotes the level of subgrid. White-dotted lines in upper and middle panels show 
the shocked regions, which indicate the first core. The elapsed time, density at the center of 
the cloud and arrow scale are denoted in each panel. 
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Fig. 2. — (a) Oblateness (e Q b) within p < 0.1p c , (b) Thick lines: Magnetic flux density (B c ; 
left axis). Thin lines: rotation period (P; right axis), (c) Magnetic field normalized by the 

1 /2 

square root of the density (B c / p c ), and (d): Angular velocity normalized by the free-fall 
timescale [f2 c /(47rGp c ) 1 / 2 ] at the center of the cloud against the central number density (n c ) 
for models SR, MR, and FR. 
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Fig. 3. — Time sequence of model SR after the protostar formation epoch . The density 
(color-scale) and velocity distribution (arrows) on the cross-section in the y = plane are 
plotted in each panel. Panels (a) through (f) are snapshots at the stages of (a) t c = 170 days 
(/ = 26), (b) 207 days (I = 26), (c) 221 days (I = 25), (d) 229 days (I = 24), (e) 236 days 
(I = 23), and (f) 265 days (I = 22), where t c denotes the elapsed time after protostar 
formation (n c = 10 21 cm -3 ). Red thick lines mean the border between the jet and the 
accretion flow (i.e., contour of v z — 0). White dotted lines in panels a-c indicate the second 
core (i.e., protostar), while white dotted line in panel / indicate the first core. The elapsed 
time after protostar formation (t c ), the elapsed time from the initial (t), gas density at the 
center of the cloud (n c ), arrow scale, and grid scale are denoted in each panel. The level of 
the subgrid is shown in the upper left corner. 
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Fig. 4. — Bird's-eye view of model SR at the same epoch as Fig. [3]f, but with different 
scales of L = 0.35 AU (left), and L = 0.087 AU (right), where L means the grid scale. The 
structure of high-density region (p > 0.1p c ; red iso-density surface), and magnetic field lines 
(black-and-white streamlines) are plotted in each panel. The structure of the jet is shown 
by purple iso-velocity surface in which the gas is outflowing from the center. The density 
contours (false color and contour lines), velocity vectors (thin arrows) on the mid-plane of 
x = 0, y — 0, and z =0 are, respectively, projected in each wall surface. 



Fig. 5. — As same as Fig. [2] but for model MR. Panels from left to right are snapshots at the 
stages of n c = 3.2 x 10 16 cm~ 3 (I = 20, panels a, b, and c), 1.2 x 10 18 cm~ 3 (I = 20; panels 
d, e, and f), and 2.3 x 10 21 cm~ 3 (/ = 20; panels g, h, and i), respectively. The white-dotted 
lines in upper and middle panels mean the first core, while the red-dotted lines in panels g 
and h indicate the protostar. 
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Fig. 6.— Time sequence of the outflow driven from the first core for model MR. The 
density (color-scale) and velocity distribution (arrows) on the cross-section in the y = 
plane are plotted in each panel. Panels (a) through (c) are snapshots at the stages of (a) 
n c = 2.5x 10 15 cm~ 3 (I = 17-19), (b) 1.5xl0 17 cm~ 3 (I = 17-19), and (c) 1.2xl0 23 cm~ 3 
(I = 17 — 19). White-dotted lines indicate the first core. Red-solid lines indicate the border 
between the outflow and accretion flow (i.e., contour of v z = 0). The elapsed time from the 
initial (t), density at the center of the cloud (n c ), arrow scale are denoted in each panel. 
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Fig. 7. — Time sequence of the jet outflow driven from the protostar for model MR. The 
density (color-scale) and velocity distribution (arrows) on the cross-section in the y = 
plane are plotted in each panel. Panels (a) through (c) are snapshots at the stages of (a) 
n c = 2.5 x 10 20 cm- 3 (I = 17-19), (b) l.lxl0 22 cm" 3 (/ = 17-19), and (c) 1.2xl0 23 cm" 3 
(I = 17 — 19). The first core and second core (or protostar) are shown by white- and red- 
dotted lines, respectively. Red-thick lines indicate the border between the jet and accretion 
flow (i.e., contour of v z = 0). The elapsed time from the initial (t), density at the center of 
the cloud (n c ), arrow scale are denoted in each panel. 
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Fig. 8. — Upper panels: the density (color-scale and white contours) and velocity distri- 
bution (arrows) on the cross-section in the y = plane for model MR at the same epoch 
as Fig. [6] c. Lower panels: the vertical component of the velocity (v z ; color-scale and con- 
tours) and magnetic field (arrows) with the same scale as each upper panel. The characters 
of 'peak' indicate the position where v z has a local peak. The first core and second core 
(or protostar) are shown by the white- and red-dotted lines, respectively. Red-thick lines 
indicate the border between the outflow and accretion flow (i.e., contour of v z — 0). 
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Fig. 9. — Bird's-eye view of model MR at the same epoch as Fig. [6]c and Fig. [7]c, but 
with different scales of L = 2.8 AU (left), and L = 1.4 AU (right). The structure of high- 
density region (p > 0.1p c ; red iso-density surface), and magnetic field lines (black-and-white 
streamlines) are plotted in each panel. The structures of the jet (v > 7kms _1 ) and outflow 
(v > 0.5 km s -1 ) are shown by purple and blue iso-velocity surfaces, respectively. The density 
contours (false color and contour lines), velocity vectors (thin arrows) on the mid-plane of 
x = 0, y = 0, and z =0 are, respectively, projected in each wall surface. 
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Fig. 10. — Time sequence of model FR. The density (color-scale) and velocity distribution 
(arrows) on the cross-section in the y = plane are plotted. Panels (a) through (c) are 
snapshots at the stages of (a) n c = 1.7 x 10 12 cm" 3 (Z = 12 — 14), (b) 2.7 x 10 14 cm -3 
(/ = 12 — 16), and (c) 5.5 x 10 14 cm -3 (Z = 16 — 17). Red-thick lines indicate the border 
between the outflow and accretion flow (i.e., contour of v z — 0). The elapsed time after the 
first core formation (t c ), elapsed time from the initial (t), density at the center of the cloud 
(n c ), and arrow scale are denoted in each panel. 



-40- 




Fig. 11. — The density (color-scale and white contours) and velocity distribution (arrows) 
on the cross-section in the y = plane (upper panels) and z = plane (lower panels) at the 
same epoch as Fig. [TOl c. but with different grid levels of I = 13 (left panels), I = 15 (middle 
panels), and I = 17 (right panels), respectively. White-dotted lines indicate the first core. 
Red-thick lines indicate the border between the outflow and accretion flow (i.e., contour of 
v z = 0). 
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Fig. 12.— Bird's-eye view of model FR (/ = 12) at the same epoch as Fig. [TOl c. The 
structure of high-density region (n > 10 12 cm -3 ; red iso-density surface), and magnetic field 
lines (black-and-white streamlines) are plotted. The structure of the outflow is shown by 
the blue iso-velocity surface inside which the gas is outflowing from the center. The density 
contours (false color and contour lines), velocity vectors (thin arrows) on the mid-plane of 
x — 0, y — 0, and z =0 are, respectively, projected in each wall surface. 
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Fig. 13.— The ratio of the toroidal to poloidal component of the magnetic field 
(i?toroidai/-B P oioidah color and contours) and velocity vector (arrows) are plotted in z > 
region for models with initially slow rotation model SR (left panels), and rapidly rotation 
FR (right panel). The plasma beta (/3 P ; color and contours) and magnetic field (arrows) are 
plotted in z < region for the same models. Red thick lines indicate the border between the 
outflow and accretion flow [i.e., the shapes of the jet (left panel) and outflow (right panel)]. 
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Fig. 14. — The ratio of each force for models SR and FR at the end of the calculation. 
Panel a: Ratio of the Lorentz to the centrifugal forces (-Piorentz /-^centrifugal) a t eacn mesh 
point in the region with v z > O.lkms -1 (i.e., the outflow and jet region). Panel b: The 
ratio of the Lorentz to the thermal pressure gradient forces (-^Lorentz/ -^thermal) in the region 
with v z > O.lkms -1 (i.e., the outflow and jet region). Panel c: The ratio of the Lorentz to 
the centrifugal forces (F Lorentz / -F centrifugal ) in the region with v z < O.lkms -1 (i.e., the inflow 
region), and Panel d: The ratio of the Lorentz to the centrifugal forces (-FL rentz/-Pcentrifugai) 
in the region with v z < O.lkms -1 (i.e., the inflow region). The horizontal axis means the 
number of the cell. The dots are the data from models FR (black dots), and SR (red dots), 
respectively. The thick lines are average of adjacent 100 cells for models FR (black lines), 
and SR (red lines), respectively. The horizontal-dotted lines indicate the value at which each 
force equals. 
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Fig. 15. — Schematic figure of the jet and outflow driven from the protostar and the fist 
core, respectively. 



